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Abstract 

We present an adiabatic quantum algorithm for the abstract problem of searching marked vertices 
in a graph, or spatial search. Given a random walk (or Markov chain) P on a graph with a set of 
unknown marked vertices, one can define a related absorbing walk P' where outgoing transitions 
from marked vertices are replaced by self-loops. We build a Hamiltonian H{s) from the interpolated 
Markov chain P{s) = {l — s)P+sP' and use it in an adiabatic quantum algorithm to drive an initial 
superposition over all vertices to a superposition over marked vertices. The adiabatic condition 
implies that for any reversible Markov chain and any set of marked vertices, the running time of 
the adiabatic algorithm is given by the square root of the classical hitting time. This algorithm 
therefore demonstrates a novel connection between the adiabatic condition and the classical notion 
of hitting time of a random walk. It also significantly extends the scope of previous quantum 
algorithms for this problem, which could only obtain a full quadratic speed-up for state-transitive 
reversible Markov chains with a unique marked vertex. 
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Introduction 

Adiabatic quantum computation was introduced by Farhi et al. in It proceeds as 
follows. Suppose that the solution of a computational problem can be encoded in the ground 
state of a problem Hamiltonian Hp. We start in the ground state of an initial Hamiltonian 
Hq which is easy to construct. Then we slowly change the Hamiltonian from Hq to Hp, so 
that the instantaneous Hamiltonian at any point in the evolution is H{s) = {1 — s)Hq + sHp, 
where < s < 1. If this is done slowly enough, then the Adiabatic Theorem of Quantum 
Mechanics guarantees that the state at all points in the evolution stays close to the 
ground state of H{s). Note that the validity of the folk version of the Adiabatic Theorem, 
as stated in many books of Quantum Mechanics such as [2], has recently been the subject 
of much debate 3|-[6|, and there is no rigorous proof of it that holds under full generality. 
It is nevertheless possible to state a more rigorous version of the theorem, so that the folk 

n 

adiabatic condition can be proved to be sufficient in many interesting cases [7|, such as the 



adiabatic version of Grover's algorithm |8l-ll0|. 



Many classical randomized algorithms rely heavily on random walks or Markov chains. 
The notion of hitting time is a useful characterization of Markov chains used when search- 
ing for a marked vertex. Quantum walks are natural generalizations of classical random 
walks. The notion of hitting time has been carried over to the quantum case in llHlTJ, 
by generalizing the classical notion in different ways. It is intimately related to the prob- 
lem of spatial search. Suppose that we are given a graph where some vertices are marked. 
Classically, a simple algorithm to find a marked vertex is to repeatedly apply some ran- 
dom walk P on the graph until one of the marked vertices is reached. The hitting time 
of P is precisely the expected number of repetitions necessary to reach a marked vertex, 
starting from the stationary distribution of P. The notions of quantum hitting time are 
based on different quantum analogues of this algorithm. They usually show some quadratic 
improvement of the quantum hitting time over the classical hitting time. However, until 
the present paper, they could only show such a relation under restricted conditions: either 
the quantum algorithm could only detect marked elements (l3|, or it could only be applied 
to state-transitive reversible Markov chains with a unique marked element 16|. Whether 
this quadratic speed-up for finding a marked element also holds for any reversible Markov 
chain and for multiple marked elements was an open question. In this paper, we answer this 
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question by the positive, by providing an adiabatic quantum algorithm for this problem. In 
addition to being more general, the algorithm is also conceptually very clean, it implements 
a simple rotation in a two dimensional subspace based on a quantum walk on the edges of 
the graph. Moreover, it reveals a close connection between the adiabatic condition and the 
notion of hitting time. 

The paper is structured as follows. In Section [T] we describe related work and in Section HTl 
we state our main result. In Section UTTl we introduce the necessary concepts such as Markov 
chains, the discriminant matrix, and the quantization of Markov chains. In Section |IV] we 
evaluate the spectrum of the interpolating Hamiltonian and in Section |V] we impose the 
adiabatic condition to calculate the running time of the adiabatic quantum algorithm. In 
Section rvT] we relate this to the hitting time of the Markov chain we started from, and show 
that the running time of the adiabatic evolution is at most the square root of the classical 
hitting time. 



RELATED WORK 



Inspired by Ambainis' quantum walk algorithm for Element Distinctness [18j, Szegedy 13 1 
introduced a powerful way of quantizing Markov chains which led to new quantum walk- 
based algorithms. He showed that for any symmetric Markov chain a quantum walk could 
detect the presence of marked vertices in at most the square root of the classical hitting time. 
However, showing that a marked vertex could also be found in the same time (as is the case 
for the classical algorithm) proved to be a very difficult task. Magniez et al. 15| extended 
Szegedy's approach to the larger class of ergodic Markov chains, and proposed a quantum 
walk-based algorithm to find a marked vertex, but its complexity may be larger than the 
square root of the classical hitting time. A typical example where their approach fails to 
provide a quadratic speed-up is the 2D grid, where their algorithm has complexity Gfn), 
whereas the classical hitting time is 0(nlogr;,). Ambainis et al. ll| and Szegedy's 13 1 



approaches yield a complexity o 



vertex. Childs and Goldstone 



19 



e 



'y/n ■ logn) in this special case, for a unique marked 
20| also obtained a similar result using a continuous-time 



quantum walk. However, whether a full quadratic speed-up was possible remained an open 



question, until Tulsi 



2l| proposed a solution involving a new technique. Magniez et al. 



extended Tulsi's technique to any reversible state-transitive Markov chain, showing that for 



3 



such chains, it is possible to find a unique marked vertex with a full quadratic speed-up over 
the classical hitting time. However, the state-transitivity is a strong symmetry condition, 
and furthermore their technique cannot deal with multiple marked vertices. It would be 
strange if one had to rely on involved techniques to solve the finding problem under such 
restricted conditions, while the classical analogue algorithm is conceptually very simple and 
works under very general conditions. 

In this paper we show that these issues can be resolved by combining the idea of the 
quantization of Markov chains with adiabatic quantum computation (note that Childs and 



Goldstone 



19 



20| showed that their algorithm for spatial search on the grid could also be 



translated into an adiabatic algorithm, but this failed to give a quadratic speed-up for low 
dimensions) . 



II. MAIN RESULT 

Before describing our quantum algorithm, let us first recall the classical algorithm on 
which it will provide a quadratic speed-up. This very simple algorithm just consists in 
walking randomly on the graph until a marked vertex is reached. More precisely, it relies 
on a Markov chain P with stationary distribution tt, and works as follows. 

Random Walk Algorithm 

1. Sample a vertex x E X according to distribution tt. 

2. Check if x is marked. 

3. If so, output X. 

4. Otherwise, update x according to P, and go back to step [21 

Let P be an ergodic Markov chain, and M be a set of marked vertices. The hitting time of 
P with respect to M, denoted by IIT(P, M), is the expected number of executions of stepH] 
during the course of the Random Walk Algorithm (where the expectation is calculated 
conditionally on the initial vertex being unmarked). Note that since the algorithm stops as 
soon as a marked element is reached, this is equivalent to using an absorbing Markov chain 
P', which acts as P on all but marked vertices, where all outgoing transitions are replaced 
by self-loops. 
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Previous attempts at providing a quantum speed-up over this classical algorithm have 
followed one of these two approaches: 

Combining a quantum version of P with a reflection through marked vertices to mimic 



a Grover operation 
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Directly applying a quantum version of P' 13|, [16 1 



The problem with these approaches is that they would only be able to find marked vertices in 
very restricted cases. We explain this by the different nature of random and quantum walks: 
while both accept a stable state, i.e., the stationary distribution for the random walk and the 
eigenstate with eigenvalue 1 for the quantum walk, the way both walks act on other states 
is dramatically different. Indeed, an ergodic random walk will converge to its stationary 
distribution from any initial distribution. This apparent robustness may be attributed to 
the inherent randomness of the walk, which will smooth out any initial perturbation. After 
many iterations of the walk, non-stationary contributions of the initial distribution will be 
damped and only the stationar y d istribution will survive (this can be attributed to the 
thermodynamical irreversibility [26| of ergodic random walks). On the other hand, this is 
not true for quantum walks, because in the absence of measurements a unitary evolution 
is deterministic (and in particular thermodynamically reversible): the contributions of the 
other eigenstates will not be damped but just oscillate with different frequencies, so that the 
overall evolution is quasi-periodic. As a consequence, while iterations of P' always lead to a 
marked vertex, it may happen that iterations of the quantization of P' will never lead to a 
state with a large overlap over marked vertices, unless the walk exhibits a strong symmetry 
(as is the case for a state-transitive walk with only one marked element, which could be 
addressed by previous approaches). 

The main new idea of our approach is that, instead of using a quantization of P or P', 
we first modify the classical random walk, and then use a quantization of the modified walk. 
The original classical algorithm consists in applying P' on the stationary distribution tt of P. 
While doing so, the system is brought far from equilibrium since vr is far (in statistical dis- 
tance) from any stationary distribution of P', which only have support on marked elements. 
The random walk will damp any non- stationary contribution of the initial distribution, but 
a quantum walk based on P or P', which is deterministic until a measurement, seems to have 
trouble with it. There is however a situation in Quantum Mechanics where contributions 



from other eigenstates will also cancel out, similarly to what happens for a random walk: 
if the system starts in a state close to the ground state of its instantaneous Hamiltonian 
H{s) {i.e., close to equilibrium), and this Hamiltonian varies slowly, the Adiabatic Theorem 
ensures that the contributions from excited states will cancel out so that the system will 
remain close to its ground state. Therefore, our strategy is to first modify the classical 
algorithm so that the system stays close to equilibrium throughout the evolution, and then 
to translate it into an adiabatic quantum algorithm. 

Consider the interpolated Markov chain P{s) = (1 — s)P + sP' (see Section [III A[) . Our 
goal is to drive the stationary distribution tt of P towards a stationary distribution of P'. 
Instead of immediately applying P' on tt, we could rather apply P{s) while slowly switching 
s from to 1, so that the system remains at all time close to the stationary distribution tt{s) 
of P{s). It can be shown that this leads to an algorithm with only a constant overhead with 
respect to Random Walk Algorithm. Therefore, this new classical algorithm still runs in 
time 0(HT(P, M)), but the difference is that at all time the system is close to equilibrium, 
so that we are in a better shape for designing a quantum analogue based on the Adiabatic 
Theorem. 

Using a Hamiltonian version of Szegedy's quantization technique, proposed by Somma 
and Ortiz [2^, we map P(s) to a Hamiltonian H{s) on the edge space Ti^Ti, where "H is the 
Hilbert space whose basis states are labeled by the vertices of the graph (see Section fill Dl) . 
The eigenstate of H{s) with eigenvalue then corresponds to the stationary distribution of 
P for s = 0, and to a distribution over marked vertices for s = 1, so that this Hamiltonian 
may be used to solve the search problem by adiabatic evolution. The algorithm consists in 
preparing the first register in the state Ivr) corresponding to the stationary distribution vr of 
P and applying the Hamiltonian H{s) using a schedule s{t) (we will specify s{t) explicitly 
later). The resulting adiabatic evolution effectively implements a rotation on the first register 
at constant angular velocity from Ivr) to a superposition over marked vertices. 

Adiabatic Search Algorithm 

1. Prepare the state |vr)|0). 

2. Apply the time-dependent Hamiltonian H{s) with sched- 
ule s{t) from t = to t = T, where T = ^v^HT(P,M). 

3. Measure the first register in the vertex basis. 
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Under the assumption that the folk adiabatic condition holds in our setup, we prove the 
following: 

Theorem 1. For any ergodic and reversible Markov chain P with a set of marked vertices 
M, the Adiabatic Search Algorithm finds a marked vertex with probability at least \ — e^ 



in a time T = ^a/HT(P, M), where HT(P, M) is the hitting time of the classical Markov 
chain P with respect to the set of marked vertices M . 

This theorem constitutes our main result and the body of this paper will be devoted 
to its proof. While it relies on the folk adiabatic condition, a similar statement can be 
made for a related quantum circuit algorithm, where no such condition is necessary, as 



shown in 



23| . Nevertheless, as explained in 23|, the intuition behind the quantum circuit 



algorithm originates from the present adiabatic quantum algorithm. 



III. PRELIMINARIES 

A. Classical interpolation 

Let us consider a Markov chain on a discrete state space X of size n, and let P be 
the n X n (row) 

chain. From now on, we will assume that the Markov chain is ergodic, meaning that it 
is both irreducible (any state can be reached from any other state by a finite number of 
steps) and aperiodic (there is no integer k > 1 that divides the length of every cycle of the 
underlying directed graph of the stochastic matrix P). Assume that a subset M C X of 
elements are marked and let m be the size of M. Let P' be the Markov chain obtained 
from P by turning all outgoing transitions from marked elements into self-loops. We call P' 
the absorbing version of P. Note that P' differs from P only in the rows corresponding to 
the marked elements (where it contains all zeros on non-diagonal elements, and ones on the 
diagonal). If we arrange the elements of X so that the marked elements are the last ones, 
matrices P and P' have the following block structure: 

P.^(PuuPum\ P':=\—-\, (1) 

\ Pmu Pmm I 
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where Puu and Pmm are square matrices of size {n — m) x {n — m) and mxm, respectively, 
while PuM and Pmu are matrices of size [n — m) x m and m x [n — m), respectively. We 
call 

P{s):={l-s)P + sP', 0<s<l, (2) 

the classical interpolation of P and P'. Note that -P(O) = P, -P(l) = -P', and has block 
structure 

- s)PmU (1 - s)PAfM + 

Moreover, note that the ergodicity of P implies that P(s) is also ergodic, except for s = 1. 



B. Stationary distribution and reversibility 

By definition, since P is ergodic, it has a unique and non- vanishing stationary distribution, 
i.e., a probability distribution vr such that ttP = vr. An ergodic Markov chain P is called 
reversible if it satisfies the so-called detailed balance condition 

Wx,y e X : -KxPxy = -n-yPyx- (4) 

This implies that for reversible Markov chains, the net flow of probability in the stationary 
distribution between every pair of states is zero. 

From now on we will assume that P is reversible. Let us argue that P{s) is also reversible. 
Let vr := (vr^ ttm) be the stationary distribution of P, where vr^/ and ttm are row vectors of 
length n — m and m, respectively. Let pm = J2xeM '^^ probability to pick a marked 

element from the stationary distribution and 7r(s) be the following distribution: 

-is) := ((1 - s)nu .m) ■ (5) 

One can check that tt{s) is a stationary distribution for P{s) for any s G [0, 1]. Moreover, for 
s G [0, 1) P(s) is ergodic and this is therefore the unique stationary distribution, while for 
s = 1 any distribution which only has support on marked vertices is stationary. Equation 
(jl]) can be used to show that 

Vs G [0,l],Vx,y G X : 7r,.(s)P,y(s) = 7r^(s)P^,(s), (6) 

which means that P(s) is reversible for s G [0, 1). Condition ([6]) is also satisfied at s = 1, 
but P(l) = P' is not ergodic, therefore stricto sensu P(l) is not reversible. 
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C. Discriminant matrix 



The discriminant matrix of a Markov chain P{s) is 

D{s) := ./Pis)oP{sy, (7) 

where the Hadamard product "o" and the square root is computed entry-wise. We prefer to 
work with D{s) rather than P{s) since a Markov chain is not necessarily symmetric while 
its discriminant matrix is. 

For a reversible Markov chain, the extended detailed balance condition implies that 
Dxy{s) = VPxy{s)Pyx{s) = Pxy{s) ^ Tixis) / Tiy{s) oi cquivalcntly 

D{s) = diag(v^) P{s) diag(v^)"\ (8) 

For s G [0, 1) the right-hand side is well-defined so that D{s) and P{s) are similar and 
therefore have the same eigenvalues. Moreover, the entry-wise square root of the stationary 
distribution ^J'K{sY is the eigenvector of D[s) with eigenvalue 1. 

At s = 1 the right-hand side of equation ([8]) is not well-defined, but it can be shown that 
both claims still hold by expanding P{s) according to equation ([3]) and considering the limit 
s — > 1. Then equation ^ becomes 

jji^l^ _ [ di^g(v^) Puu diag(v/7iv) ^ 

This implies that D{1) is similar to P := ( ? ) , and in turn to P(l) = ( ) 

To see that \Jt{{1Y is an eigenvector of D{1) with eigenvalue 1, note that 7r(l) = (0 'n'M)/pM, 

and -D(l) acts as the identity on marked elements (this follows from equations ([5]) and (IHl), 

respectively). 




D. The quantum Hamiltonian 

Szegedy [l^ proposed a general method to map a random walk to a unitary operator 
that defines a quantum walk. Recently Somma and Ortiz 22| showed how Szegedy's method 
may be adapted to build a Hamiltonian. We apply this method to the random walk P{s). 
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The first step is to map the rows of P{s) to quantum states. Let "H be a Hilbert space of 
dimension n = \X\. For every x G X we define the following state in T-L: 



yex 



Following Szegedy 



13| . we then define a unitary operator V{s) acting on ® H as 



V{s)\x,0) := \x)\p,{s)) = J2 \IP.y{s)\x,y), (11) 

when the second register is in some reference state |0) G "H, and arbitrarily otherwise. 

Let S be the gate that swaps both registers. When restricted to |0) in the second register, 
the operator V'^{s)SV{s) acts as D{s): 

{xMV\s)SV{s)\y,Q) = {x\py{s)){p^{s)\y) = sJrJs)P^) = D^y{s). (12) 



Following 22|], we now define the Hamiltonian H{s) on "H "H as 

H{s) ■.= i[V\s)SV{s),n,], (13) 

where Ho := /(8)|0)(0| is the projector that keeps only the component containing the reference 
state |0) in the second register and [A, B] := AB — BA is the commutator. 

IV. SPECTRAL DECOMPOSITION OF H{s) 

To understand the properties of the Hamiltonian -ff(s), let us find its spectral decompo- 
sition. We will relate its spectrum to that of D{s). 

A. Diagonalization of D{s) 

Recall from equation ([7j) that D{s) is real and symmetric. Therefore, its eigenvalues are 
real and its eigenvectors form an orthonormal basis of l-L with real amplitudes. Let 

D{s) = Y,Hs)\v,{s)){v,{s)\ (14) 

i=l 

be the spectral decomposition of D{s). We can make the eigenvalues of P{s) and hence also 
of -D(s) to be non-negative by replacing P with (P + /)/2. Note that this will only modify 
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the hitting time of the Markov chain by a factor of 2. Hence, from now on without loss 
of generahty we assume that all eigenvalues of D{s) are non- negative. In addition, we can 
arrange them so that 

< Ai(s) < A2(s) < ■■■ < A„(s). (15) 

From the Perron-Frobenius theorem we have that Vz : Aj(s) < 1 and A„(s) = 1. In 
addition, for any s G [0, 1) the Markov chain P{s) is ergodic and "^i ^ n : Aj(s) < 1. On the 
other hand, at s = 1 the Markov chain is not ergodic and has eigenvalue 1 with multiplicity 
m. We may summarize this as follows: 

A„_i(s) < A„(s) = 1, VsG[0,l), (16) 
A„_^(l) < A„_„+i(l) = ■ ■ ■ = A^„(l) = 1. (17) 

Recall from Section UlI CI that a/vt^s)^ is an eigenvector of D(s) with eigenvalue 1 for any 
s G [0, 1], so let us choose \vn{s)) in equation ( |T^ as 

\vn{s)) := (18) 



B. Diagonalization of H{s) 

Now, let us express the eigenvalues and eigenvectors of H{s) in terms of those of D{s). 
First, let us break up the Hilbert space 'H^'H into mutually orthogonal subspaces that are 
invariant under H[s). Let 

. Bk{s) := span{\vk{s),0),V^{s)SV{s)\vk{s),0)}, Vfc ^ n, 
• Bn{s) := span{|f„(s),0)}. 

Note that {vkis),0\V\s)SV{s)\vkis),0) = {vkis)\D{s)\vkis)) = Afe(s) by equations and 
( |T4l) . Recall that for s < 1, we have Xk{s) ^ 1 for any k ^ n and thus 

V\s)SVis)\vk{s),0) = Xkis)\vk{s),0) + ^l-h{sy\vk{s),0)^ (19) 



for some [28| unit vector \vk{s), 0)-*- orthogonal to \vk{s), 0) and lying in the subspace Bk{s). 
We also define by continuity |f a;(1), 0)-*- := lims_!>i \vk{s),0)-^. 

Following Somma and Ortiz, who were in turn relying on Szegedy's work, we may now 
characterize the spectrum of H[s). 
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Lemma 1 ( 13|, |22|). H{s) accepts the following eigenvalues and eigenstates. 



• On Bk{s), Wk 7^ n: 

EkKs) ■.= ±^l-\k[s)\ |^^(s)):= , (20) 

• On Bn{s): 

Enis):=0, l^nis)) ■.= \Vnis),0). (21) 

• OnB^is): 

F,(s):=0, \^,{s)), je{l,...,(n-l)2}, (22) 
where {|$j(s))} defines an arbitrary basis of B^{s). 

Proof. We consider the case s G [0, 1), the case s = 1 follows by continuity. Since 
V''{s)SV{s)\vn{s),0) = D{s)\vn{s))\0) = \vn{s),0), wc immediately have that |f„(s),0) is 
an eigenstate of H{s) with eigenvalue 0. For k ^ n, note that 

V\s)SV{s)Uo\vk{s),0) = Xk{s)\vk{s),0) + ^1 - Xk{sy\vk{s),0)^ , (23) 
UoV\s)SV{s)\vk{s),0) = Xk{s)\vk{s),0}. (24) 

By combining these expressions we get 

His)\vkis),0) = t^/l-Xkis)^\vkis),0)^, (25) 
H{s)\vj,{s),0)^ = -t^/l-Xk{sy\vk{s),0), (26) 

where the second line follows from the fact that H{s) is Hermitian and traceless. In other 
words, H{s) acts on subspace Bk{s) as a/1 — Xk{s)'^ay, where ay := is the Pauh y 

matrix, which yields equation (!20|) . 

Since B^{s) is the orthogonal complement of the union of invariant subspaces, it is also 
an invariant subspace for H{s). Note that H{s) restricted to this subspace is equal to zero, 
hence the remaining — (2n — 1) = (n — 1)^ eigenvalues of H{s) are zero. □ 



V. THE QUANTUM ADIABATIC THEOREM 



In adiabatic quantum computing it is a common practice to associate the intermediate 
state of the computation with the ground state {i.e., the lowest energy eigenstate) of the 
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Hamiltonian. However, in our case the spectrum of H{s) is symmetric about zero and the 
state that we are interested in hes in the middle of the spectrum. Hence, we will not use the 
ground state of H{s), which has negative energy, but instead we will use the zero-eigenvector 
given in equation fl^ . Indeed, equation flTS]) shows that this state is closely related 
to the stationary distribution tt{s) of P{s). In particular, the problem would be solved if 
we can reach the state |\1/„(1)), as measuring the first register of this state yields a vertex 
distributed according to 7r(l), which only has support on marked vertices. 

A. The adiabatic condition 

We initially prepare the system in the zero-eigenvector |\l/„(0)) of H{0) and then start to 
change the Hamiltonian H{s) by slowly increasing the parameter s from to 1 according 
to some schedule s{t). If the schedule s{t) is chosen so that it satisfies certain conditions, 
the system is guaranteed to stay close to the intermediate zero-eigenstate Then, at 

s = 1, the state will be close to = |t'n(l), 0), where the first register only has overlap 

over marked vertices, so that a measurement yields a marked vertex with high probability. 
Note that in our case the zero-eigenspace -B„(s) UB-^{s) of the Hamiltonian H{s) has a huge 
dimension, so we have to make sure that the non-relevant part B^{s) is totally decoupled 
from (the only zero-eigenvector that is relevant for our algorithm) before we apply 

the adiabatic condition. In particular, we want to show that 

{H-^)\-j^\^nis)) = (27) 

for any j G {1, . . . ,{n — 1)^}, since this would imply that during the evolution |\I^„(s)) is 
not leaked into the subspace B^{s) spanned by states |$j(s)). To see that this is indeed the 
case, note that 

|$,(s)) L span {|vl/±(,)), . . . , |vl/±_,(s)), {^^(s))} (28) 

for any j G {1, . . . , (n — 1)^}, since the eigenvectors of H{s) form an orthonormal basis. In 
particular, 

|$,(.)) ± span {|M/+(s)) + |vl/-(s)), . . . , + l^^-il^)), ^8))} 

= span {|t;i(s), 0), . . . , |t;„_i(s), 0), |t;„(s), 0)} . (29) 
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Recall from equation fl2Tl) that ^|\E'„(s)) = ^|^^n(s))|0). Hence, the inner product in equation 
(!27|) indeed vanishes. Thus, we can safely apply the adiabatic condition only for the relevant 
subspace in which the zero-eigenstate is not degenerate. 

The folk version of the Adiabatic Theorem [2] states that during the evolution the state 
of the system li'it)) is guaranteed to stay close to the intermediate zero-eigenstate |\E'„(s)), 
more precisely, 

Vt: \{^n{s{t))\^P{t))\'>l-e^, (30) 
as long as the adiabatic condition 



{EM - EJs))' 



is satisfied. While this condition is known not to be sufficient in full generality (see e.g. the 



discussion in j?!), we will assume that it can be applied in our setup. We will discuss about 
how this assumption may be suppressed in Section IVIIi 

If we insert all eigenvalues and eigenvectors from equations (1201) and fl2T|) . the adiabatic 
condition (1311) can be written purely in terms of the eigenvalues and eigenvectors of the 
discriminant matrix D{s): 

W: E l<"-W^f;j'»l' <.-. (32) 
B. Rotation in a two-dimensional subspace 

In this section we will provide a simple interpretation of the evolution of the eigenvector 
\vn{s)). First, let us define the following superpositions over all elements, marked elements, 
and unmarked elements, respectively: 

k):=5Zv^l^)' l^):=^$^v^k), 1^):= rr^ : ^ v^|x). (33) 

Now, we show that \vn{s)) is subject to a rotation in the two-dimensional subspace TZ := 
span{|[/), |M)}. From equations ([5]) and ( |T8|) . we obtain 

l^n(s)) = ^/7Tx{s)\x) = ^ == I Vl - S ^ y/^\x) + ^ y/^\x) ) . (34) 
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Using the definition of \U) and \M) we can rewrite tliis simply as 

\v^[s)) = cose{s)\U) + sme{s)\M), (35) 

wliere 



^(s):=arcsin^^— (36) 

Let us clioose tlie scliedule s{t) so tliat tlie evolution of \vn{s)) as defined by equation 
( 135|) corresponds to a rotation witli constant angular velocity in the subspace TZ. If T is the 
length of the evolution and s : [0, T] — )■ [0, 1] is defined as 

s(t) := - — ^ ( 1 . , n \ I ' ^0 := arcsin y/p^, ■= ^ arccos ^/p]^, (37) 

then 

e{s{t)) = ujt + eo. (38) 

Let us choose a vector |f^(s)) such that {|f(s)), |f;^(s))} is an orthonormal basis of TZ 
for every s: 

\v:^{s)) := - sin e{s)\U) + cos e{s)\M). (39) 
Then from equations ( |35|) and (!38|) we easily find that 

^l^n(s(t))) = u\v^{s{t))) = ^ arccos ,/jJ^\v^{s{t))). (40) 

Note that arccos ^/pm < f • Therefore, we can rewrite the adiabatic condition ( 132|) as follows: 



""''■aT^ ^ i-xlis) 

If this condition is satisfied, equation ( 130|) implies that we obtain at time t = T a state 
|'?/'(7')) close to |\E'„(1)) = |t;„(l))|0) = |M)|0), so that measuring the first register yields a 
marked vertex with probability at least 1 — e^. 

VI. RUNNING TIME OF THE QUANTUM ALGORITHM 
A. Choice of running time T 



We have to change the parameter s slowly for the evolution to be adiabatic, thus we want 
to choose T big enough so that condition (j4T|) holds. Recall from Section IIV Al that we can 
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assume that Xk{s) > 0. Thus, 1 — A^(s) = (l + Afc(s)) (l — \k{s)) > 1 — Xk{s). Let us impose 
a shghtly stronger condition on T in equation (HTl) by replacing 1 — A^(s) with 1 — Afc(s). 
In addition, let us choose the smallest T that still satisfies the inequality and use it as the 
running time of our adiabatic algorithm: 



71 

T := — max 

2e o<s<i 



V 1-Xk{s) ■ ^ ^ 



It turns out that there is an interesting relationship between this quantity and the hitting 
time of the Markov chain P. 



B. Hitting time of a Markov chain 

Let us first give an explicit expression for the hitting time HT(P, M) as defined in Sec- 
tion Ull Let Om and \m (resp., 0[/ and \u) be the all-zero and all-one row vectors of dimension 
m (resp., n — m). Furthermore, let t^m '■= t^m/Vm and VTf/ := 71(7/(1 — pm) be the row vec- 
tors describing distributions over marked and unmarked vertices. Then, the distribution of 
vertices at the the first execution of step |4] of Random Walk Algorithm is {liu Oa/), and 
from the definition of the discriminant -D(l), we have 

oo 

HT(P, M) := ^^Pr[No marked vertex found after t steps from {txu Om)] (43) 
t=o 

oo 

= ^m)P\1){Iu ^mY (44) 

t=0 

oo 

= Y.{U\D\l)\U), (45) 

t=o 

where the last equality follows from equation ([9]). We will show that the running time of 
our adiabatic quantum algorithm is directly related to the square root of the hitting time 
HT(P, M). In order to do this, we define the following extension of the hitting time. Let 

oo 

HT(.) := Y.^U\{D\s) - \u^{s))Ms)\)\U). (46) 

Note that HT(1) = HT(P, M) since ([/|f„(l)) = 0. This justifies to consider HT(s) as an 
extension of the hitting time. Intuitively, HT(s) may be understood as the time it takes for 
P(s) to converge to its stationary distribution, starting from (tt^/ Om)- For s = 1, the walk 
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P(l) = P' converges to the (non-unique) stationary distribution (Of/ ttm), which only has 
support over marked elements. 

Using the expansion (1 — x)~^ = ^^g^* spectral decomposition (fUj) of the 

discriminant D{s), we have 

oo 

HT(5) = EE^l-(^)(f^l^^(^))(^'^-(^)|f^) (47) 

fc^n t=0 

~ ^ 1-Ms) • ^^^^ 

C. Relation between the running time and the extended hitting times 

Let us express the running time T from equation (H2|) in terms of HT(s). Define 

Note that both T and HT(s) can be expressed in terms of A{s) as follows: 

T = ^ max ^{v^{s)\A{s)\v^{s)), HT(.) = {U\A{s)\U). (50) 

By definition, we have A{s)\vn{s)) = 0, which together with equation ( 135|) implies that 

^(^)|M) = -^^A(s)|f/). (51) 

Using this and the definition of \v^{s)) in equation (l39l) . we see that {v:^{s)\A{s)\v^{s)) = 
{U\A{s)\U) / sin^ 9{s). Thus we get the following relationship between T and HT(s): 



T = - max (52) 

2e o<s<i sin 61(5) 

To relate T and the usual hitting time HT (P,M) of P, we first provide an explicit 
expression for HT(s) in terms of HT(P, M) (the proof is given in the appendix). 

Lemma 2. HT(s) = HT(P,M) sin^^(s). 

Now, recalling the definition of 6{s) in equation (!36|) . it is easy to see that the maximum 
in equation fl52l) is attained at s = 1. This immediately implies that the running time of the 
adiabatic quantum algorithm is given by 

TT 



T = -v/HT(P,M), (53) 

therefore providing a quadratic improvement over the classical hitting time. This also con- 
cludes the proof of Theorem [TJ 
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VII. CONCLUSION AND DISCUSSION 



Our quantum algorithm defines a new notion of quantum flitting time, which is quadrat- 
icaify smaffer than the cfassicaf hitting time for any reversible Markov chain and any set of 
marked elements. While previous approaches were subject to various restrictions, e.g., the 
quantum algorithm could only detect the presence of marked elements jl^, did not always 



provide full quadratic speed-up 



15|, or could only be applied for state-transitive Markov 



chains with a unique marked element 16|, our adiabatic approach only requires minimal 
assumptions. Indeed, it can be shown that the only remaining condition, reversibility, is 
necessary. Let us consider the Markov chain on a cycle P = ^-y^, where C implements a 
clockwise shift, i.e., C\x) = |(x + 1) mod n). This Markov chain is ergodic but not re- 
versible. While its classical hitting time is of order Q{n), a simple locality argument implies 
that any quantum operator acting locally on the cycle requires a time Q{n) to find a marked 
vertex, so that a quadratic speed-up cannot be achieved. Magniez et al. 16|] have also 
shown that under reasonable conditions the quadratic speed-up is optimal. This provides 
evidence that our result is both as strong and as general as possible. 

While our result relies on the assumption that the folk adiabatic condition is sufficient, 
this assumption could be suppressed in different ways. One option would be to actually 

prove that the folk Adiabatic Theorem holds in our setup, as was previously done for the 

i 

adiabatic version of Grover's algorithm \J\. Another option would be to circumvent adiabatic 
evolution altogether, by using the phase randomization technique of Boixo et al. [24] . Their 
technique provides a quantum circuit realizing the same evolution as the adiabatic approach 



with a similar running time, but without re 
the quantum circuit algorithm described in 



ying on the adiabatic condition. This leads to 



23|. 



Finally, note that in order to design the schedule our algorithm requires to know Pm 
and the order of magnitude of HT(P, M). These assumptions are standard in other quantum 
algorithms for this problem. In particular, a similar issue arises in Grover's algorithm when 
the number of marked elements is unknown. In Grover's case, there are techniques to deal 
with this issue [25|, and similar techniques could be apphed in our case. While we do not 
provide a full answer to these questions in the present paper, they do not present any new 
technical difficulty and we refer the reader to 23| where a full study of these implementation 
issues is provided for a related quantum circuit algorithm. 
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Appendix: Proof of Lemma [2] 

In this section, we will often use /(s) as a shorthand form of ^/(s). We will show that 
the derivative of HT(s) satisfies the following lemma. 

Lemma 3. The derivative o/HT(s) is related to HT(s) as: 

Note that Lemma [2] follows directly from there, since HT(s) = HT(P, M) sin^ 6{s) sat- 
isfies the differential equation for HT(s) with the boundary condition HT(1) = HT(P, M). 
Therefore, it remains to prove Lemma [3l 

Before proving Lemma [3], let us consider the derivative of the discriminant D{s). Let 
= J2x£M t>e the projector onto the m-dimensional subspace spanned by marked 

vertices, and let {X, Y} := XY + YX be the anticommutator of X and Y. 
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Lemma 4. D{s) = j^^{Um, I - D{s)} . 

Proof. Note from equation ([3]) that D{s) has the following block structure: 



v^w^ v(i-.')(^™op;„) ^^^^ 



Hence, 



^(^)=| , ^ ^-^^-^'^"^ 7^1- (A.3) 



Observe that D{s) + D{s)} = -^Um, which implies the lemma. □ 

We are now ready to prove Lemma [3l 

Proof of Lemma 0. In this proof we will often omit to write the dependence on s explicitly. 
From equation (1501) we have 

^RT = {U\A\U). (A.4) 

Note that A = B'^ -U^, where 

P:=/-D + n„, U„.:=K){Vr,\. (A.5) 

For any invertible matrix M we have ^(M~^) = — M~^MM~^. Therefore, 

A = -B-\-D + tln)B-^^tln. (A.6) 

Hence, we have ^ HT = hi + h2 + h^, where 

h := {U\B-^DB-^\U), (A.7) 
h2:=-{U\B-'Ur,B-'\U), (A.8) 
h := -{U\Un\U). (A.9) 

Let us evaluate each of these terms separately. We can use Lemma H] and the definition of 
B to express the first term hi as follows: 

2(1 = {U\B-^{Um,I - D}B-^\U) (A.IO) 

= {U\B-'{{IIm,B} - {Um,I14)B-'\U) (A.ll) 
= {U\{B-\Um}\U) - {U\B-'{Um,II,,}B-'\U). (A.12) 
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Note that IIm\U) = 0, thus the first term vanishes. Also note that B ^\vn) = \vn) and B ^ 
is Hermitian, thus 

2(1 - s)hi = -2{U\B-^UMlln\U). (A.13) 

Note from equation fl5S]) that UmIvu) = sin6'|M). Moreover, using equation fl^ we can 
simphfy equation (IA.13P even further: 

2(1 -s)/ii = -2cosesme{U\B-^\M) (A.14) 

= -2cos^sin^(f/|(A + n„)|M) (A. 15) 

= -2 cos e sme{U\A\M) - 2 cos^ 9 sin^ 9 (A.16) 

= 2cos2^((f/|A|t/) -sin^^). (A.17) 

Let us now consider the second term 

h2 = -{U\B-^UnB-^\U) (A.18) 

= -cose{{v^\B-'\U) + {U\B-'\vn)) (A.19) 

= -2cose{U\{A + Un)\vn) (A.20) 

= -2^cos^(-sin0(t/|A|f/) +cos^(t/|A|M)) (A.21) 

= -2^cose (-sm0(U\A\U) - ^^^^(U\A\U)] (A.22) 
\ sin / 

= 2e--{U\A\U), (A.23) 
sm t7 

where we have used the fact that 

\u^) = e{-sme\u) + cose\M)), (A.24) 

and Ilnl-Un) = 0. Similarly, for the last term we have 

/i3 = -(f/|Il„|f/) = 2^cos^sin^. (A.25) 

Putting all the terms back together, we have 

{U\A\U) = ({U\A\U)-sm^e) + 26— -{U\A\U) + 26 cos 6 sin 6 (A.26) 

1 — s sin 6 

cos6 ( ■ cos6'sin6'\ ^ ^ / a cos6'sin6'\ 
26 + — l^V\A\V) + cos 6 s\Ty6 26 . 



sin 6 \ \ — s ^ 

From the definition of 6 in equation f l36|) it is straightforward to check that 

• cos 6* sin 6* / . x 

26 = — — — (A.27) 

which together with equations (]A.4I) and (IA.26I) implies the lemma. □ 
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